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ABSTRACT 

We report a detection of the excited 2 20 — 2 n rotational transition of para-H 2 in APM 08279+5255 
using the IRAM Plateau dc Bure interferometer. At z = 3.91, this is the highest-redshift detection 
of interstellar water to date. From LVG modeling, we conclude that this transition is predominantly 
radiativcly pumped and on its own does not provide a good estimate of the water abundance. However, 
additional water transitions are predicted to be detectable in this source, which would lead to an 
improved excitation model. We also present a sensitive upper limit for the HF J = 1 — absorption 
toward APM 08279+5255. While the face-on geometry of this source is not favorable for absorption 
studies, the lack of HF absorption is still puzzling and may be indicative of a lower fluorine abundance 
at z = 3.91 compared with the Galactic ISM. 

Subject headings: cosmology: observations - galaxies: active - galaxies: high-redshift - galaxies: 
individual (APM 08279+5255) - ISM: molecules 



1. INTRODUCTION 

Molecules such as CO or HCN have been commonly 
used as tracers of molecular gas in high-redshift galaxies. 
However, rec ent observations with the Herschel Space 
Observatory (|Pilbratt et al.l 120101) have shown strong 
spectroscopic signatures from other light hydrides, such 
as water, H2O" 1 ", or HF, in nearby active galaxies (e.g., 
van der Werf et"afll20iq ). These lines are blocked by the 
Earth's atmosphere, but can be observed, redshifted, in 
distant galaxies using the current millimeter and sub - 
millimeter facilities. For example, lOmont et al.l (|2011[ ) 
have recently reported a detection of water in J090302- 
014127B (SDP.17b) at z = 2.30. 

One of the excitin g recent results from HIFI 
(|de Graauw et al1l2010f ) is the detection of widespread 
absorption in the fundamental J = 1 — rota- 
tional transi tion of hydro g en fluoride toward Galac- 
tic sources (INeufeld et al l l2010t IPhillips et al.l l2010t 
ISonnentrucker et al.l 120101: iMonie et al.l 120111 ). Fluorine 
is the only atom that reacts exothermica lly with H 2 
(jNeufeld et al.ll2005tlNeuTeld fc Wolfirel2"009h . The prod- 
uct of this reaction, HF, is thus easily formed in re- 
gions where H 2 is present and its very strong chemical 
bond makes this molecule relatively insensitive to UV 
photodissociation. As a result, HF is the main reser- 
voir of fluorine in the interstellar medium (ISM), with a 
fractional abundance of ~3.5 x 10 -8 relative to H 2 typ- 
ically m easured in diffuse molecular clouds within the 
Galaxy (INeufeld et al.l l2010t ISonnentrucker et al.1 120101 : 
IMonie et al.ll201lD . 
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Inte rstellar HF was first detected by INeufeld et al.l 
l|1997f) with the Infrared Space Observatory (ISO). The 
,7 = 2 — 1 rotational transition was observed in absorption 
toward Sagittarius B2, at a low spectral resolution using 
the Long- Wavelength Spectrometer (LWS). The HIFI in- 
strument allows for the first time observations of the fun- 
damental rotational transition of HF at 1.232476 THz 
to be carried out, at high spectral resolution. Given the 
very large Einstein A coefficient (2.423 x 10 -2 s _1 ; critical 
density ^3 x 10 9 cm -3 ), this transition is generally ob- 
served in absorption against dust continuum background. 
Only extremely dense regions with strong IR radiation 
field could possibly generate enough collisional or radia- 
tive excitation to yield an HF feature with a positive 
frequency-integrated flux0 The HIFI observations cor- 
roborate the theoretical prediction that HF will be the 
dominant reservoir of interstellar fluorine under a wide 
range of interstellar conditions. The HF J = 1 — transi- 
tion promises to be a excellent probe of the kinematics of, 
and depletion within, absorbing material along the line 
of sight toward bright continuum sources, and one that is 
uncomplicated by the collisionally-excited line emission 
that is usually pres ent in the spectra of other gas trac- 
ers. As suggested bv lNeufeld et al.l (|2005| ) . redshifted HF 
,7=1 — absorption may thus prove to be an excellent 
tracer of the interstellar medium in the high-redshift Uni- 
verse, although only the gas reservoir in front of a bright 
continuum background can be studied by means of the 
HF absorption spectroscopy. 

Water is another interstellar molecule of key impor- 
tance in astrophysical environments, being strongly de- 
pleted on dust grains in cold gas, but abundant in warm 
regions influenced by energetic process assoc iated with 
star formation (see (van Dishoeck et all 120111 and refer- 
ences therein). The excited 2 2 o — 2n transition of p-H 2 0, 
with a lower level energy of 137 K, has a frequency of 



5 HF emission has recently been reported in the ex- 
treme environment of the nearby ultraluminous galaxy Markar- 
ian 231 in Herschel/SPIR E observations at low spectral resolution 
II van der Werf et al.|[20TOl'l ■ 
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1.228788 THz and can be observed simultaneously with 
the J = 1 — transition of HF in high-redshift systems. 
Consequently, we have searched for the HF J = 1 — and 
H2O 220 — 2n transitions, rcdshiftcd down to 251 GHz, 
in APM 082791+5255 using the IRAM Plateau de Bure 
interferometer. 

The broad absorption line (BAL) quasar 
APM 082791+5255 at 2=3.9118, with a true bolo- 
metric luminosity of (0.7 — 3) x 10 14 L Q . is one of the 
most luminous objects in the Universe (jDownes et al.l 
I1999T ) . CO lines uptoJ = ll — 10 have been detected 
using the IRAM 30-m telescope. IRAM PdBI high 
spatial resolution observations of the CO J — 4 — 3 and 
9 — 8 lines, a nd of the 1.4 mm d ust continuum have been 
presented by IWeifS et al.l (120071 ) . The line fluxes in the 
CO ladder and the dust continuum fluxes are well fit by 
a two-component model that invokes a cold component 
at 65 K with a high density of n(H 2 ) = 1 X 10 5 cm -3 , 
and a warm, ~ 220 K, component with a density of 
1 x 10 4 cm 3 . Weifi et al. argue that the molecular lines 
and the dust continuum emission arise from a very com- 
pact (r ~ 100 — 300 pc), highly gravitationally magnified 
(m ~ 60 — 110) region surrounding the central AGN. 
Part of the difference relative to other high-z QSOs may 
therefore be due to the configuration of the gravitational 
lens, which gives us a high-magnification zoom right into 
the central 200-pc radius of APM 08279+5255 where IR 
pumping plays a significant role for the excitation of the 
molecular lines. 

High-angular resolution (0.3") VLA observations of 
the CO J = 1- emission in APM 08297+5255 
(|Riechers et al.ll2009[) reveal that the molecular emission 
originates in two compact peaks separated by <0"4 and 

is virtually co-spatial with the optical/near infrared con- 
tinuum emission of the central active galactic nucleus 
(AGN). This morphological similarity again indicates 
that the molecul ar gas is located in a compact region, 
close to the AGN. IRiechers et al.l (|2009l ) present a revised 
gravitational lens model of APM 08297+5255, which in- 
dicates a magnification by only a factor of 4, in contrast 
to much higher magnification factors of ~100 suggested 
in earlier studies. Their model suggests that the CO 
emission originates from a ~550 pc radius circumnuclear 
disk viewed at an inclination angle of <25°, or nearly 

f ace-on. The to tal m olecular mass is then 1.3 x 10 11 M Q . 

iWeiB et al.l (|2007f ) first pointed out the impor- 
tance of infrared pumping for the excitation of HCN 
in APM 08279+52 55. Subsequent observations of 
IRiechers et al.l (|2010D reveal surprisingly strong J = 6 — 5 
emission of HCN, HNC, and HCO + in the host galaxy, 
providing additional evide nce that these t r ansiti ons are 
not collisionally excited. IRiechers et al.l (|2010l) argue 
that the high rotational lines of HCN can be explained 
by infrared pumping at moderate opacities in a ^220 K 
warm gas and dust component. These findings are con- 
sistent with the overall picture in which the bulk of the 
gas and dust is situated in a compact, nuclear star burst, 
where both the AGN and star formation contribute to 
the heating. 

Prior to the observations reported here, water had 
not been detected in APM 08279+5255. However, 
IWagg et all (f2006h give an upper limit of 0.7 Jykms 
(3tr) for the ground state lio — Iqi ortho-H 2 line. 



2. OBSERVATIONS 

Observations of APM 08279+5255 presented here were 
carried out on 2010 June 22, September 21-22, and De- 
cember 15, using the Plateau de Bure interferometer. 
Visibilities were obtained in the CD set of configurations 
of the six-element array in June and December and with 
a four-element subarray in September, totalling 4.9 hr of 
on-source observations. 

Data reduction and calibration were carried out using 
the GILDAS software package in the standard antenna 
based mode. The passband calibration was measured 
on 3C454.3, and amplitude and phase calibration were 
made on 0749+540, 0836+716 and 0917+449. The ab- 
solute flux calibration, performed using MWC349 as the 
primary calibrator (2.55 Jy at 250 GHz), is accurate to 
within 10%. Point source sensitivities of 4.5 mJy beam -1 
were obtained in channels of 20 MHz, consistent with the 
measured system temperatures (200-300 K). The con- 
version factor from flux density to brightness tempera- 
ture in the 1.6" x 1.3" (PA=22°) synthesized beam is 
9.1 K(Jybeam -1 ) -1 . 

3. RESULTS 

Figure [T] (upper panel) shows a spectrum of 
APM 08279+5255 near the rest frame frequency of the 
HF J = 1 — transition, integrated over the PdBI image, 
which also covers frequencies of the para-H20 2 2 o — 2xi 
and HCN J = 14- 13 transitions, in addition to HF. The 
continuum is detected with a high SNR. The integrated 
flux density, computed from emission free channels, is 
34+0.55 mJy, consi stent with the pre vious measurements 
of the source SED (jWeifi et alJl2007t ). No HF absorption 
is seen, with a 3<r upper limit of 1.5 Jykms -1 , assuming 
a FWHM line width of 500 kms -1 , as implied by earlier 
CO observations. 

The integrated line and continuum fluxes given above 
impose a 3cr upper limit of 0.092 for the velocity-averaged 
HF J = 1 — optical depth (velocity-integrated optical 
depth rdv = 46 kms -1 ). The corresponding column den- 
sity of cold HF in front of the continuum source ca n then 
be computed using eq. (3) of iNeufeld et al.l ([2010;) to be 
AT(HF) < 1.1 x 10 14 cm -2 . Given the typical Galactic 
HF/H2 abundance ratio of ~3.5 x 10 -8 , this value would 
imply an average H2 column density ^3x 10 21 cm -2 lying 
in front of the continuum source in APM 08279+5255. 
This value is three orders of magnitude below the beam 
averaged H2 column density inferred from the dust con- 
tinuum flux observed toward the source^ 

The para-H20 220 — 2n line is clearly detected with 
the integrated line flux density of 6.7 + 0.8 Jykms -1 . A 
Gaussian fit gives a line width of 510 ± 70 kms -1 , con- 
sistent with that of CO. Figure [1] (lower panel) shows 
spatial distribution of the para-H20 220 — 2n emis- 
sion (red contours) superposed on a grayscale image of 
the dust continuum. The line a nd continuum emissio n 
peak toward sources A and B of IRiechers elaTI (f2009h . 

6 I n their two-compo nent model for the dust continuum radi- 
ation, IWeifi et afl 120071 derived estimates of ~2.8 X 10 9 m -1 M 
and 680 pc, respectively, for the dust mass and magnified radius 
for APM 08279+5255, where m is the lens magnification. These 
values imply an average H2 column density of 1.2 X 10 25 cm -2 for 
an assumed dust-to-gas mass ratio of 100. Approximately 50% of 
this gas should be in front of the continuum source. 
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Fig. 1.— (Top) Integrated spectrum of APM 08279+5255. Ver- 
tical arrows mark redshiftcd frequencies of para-H20 220 — 2n, 
HF J = 1 — 0, and HCN J = 14 - 13, assuming z=3.9118. Veloc- 
ity scale is with respect to the HF J = 1 — frequency. HF and 
HCN lines are not detected. (Bottom) Distribution of the velocity- 
integrated para-H 2 2 2 rj - 2n line intensity in APM 08279+5255 
(red contours) superposed on a grayscale image of the 1.23 THz 
(rest frame) dust continuum emission. Contour levels are -2, 2, 
3, 4, 5, 6, 7 times the rms of 0.8 Jykms -1 beam -1 . White sym - 
bols mark the locations of sources A and B (Ricchcrs ct al. 2009). 
Synthesized beam is show as a white ellipse in the lower-left corner. 

The small offset between the H 2 and continuum emis- 
sion is not significant at the spatial resolution of the 
present observations. Implication for water excitation 
in APM 08279+5255 are discussed below. 

Some excess emission above the continuum level is seen 
near the frequency of the HCN J = 14—13 line, how- 
ever, the result does not constitute a detection at the 
sensitivity limit of the present observations. 

4. WATER EXCITATION 

In modeling the water line flux observed from 
APM 08279+5255, we have computed the CO and H 2 
line luminosities expected for an isothermal, constant 
density medium. We solved the equations of statistical 
equilibrium for the H 2 and CO level populations, mak- 
ing the large velocity gradient (LVG) approximation and 
treating the effects of radiative trapping with an escape 
probability metho d. We adopted the rate coeffi cients of 
lYang et all (|2010l) and iFaure fc Josselinl (|2008l) . respec- 
tively, for the excitation of CO and H2O in collisions 
with H2, and we assumed an ortho-to-para ratio (OPR) 



of 3 for both H 2 and H 2 0. Following IWeifi et all (12007ft . 
we neglect any effects of dust extinction upon the emer- 
gent CO and H 2 line fluxes; although the dust optical 
depths at THz frequencies raise the possibility that such 
effects could be significant, their importance depends 
strongly on the geometry of the source and the spatial 
relationship between the warm dust and the molecular 
emission region. In Galactic hot cores with H 2 column 
densities comparable to that in APM 08279+5255, the 
para-H 2 2 2 o — 2n line can be seen with net-emission 
flux (Orion KL, NGC6334I), or in absorption (Sagittar- 
ius B2), depending on the specific source geometry. A 
mixture of such regions may contribute to the observed 
spectrum of APM 08279+5525, leading to partial can- 
cellation of the emergent line flux. Using the relative 
strengths of the multiple CO transitions observed by 
IWeifi et all (|2007D to constrain the gas temperature, den- 
sity, and velocity gradient, we thereby obtained as best fit 
parameters the values T = 10 2 ' 35 K, n(H 2 ) = 10 42 cm" 3 , 
and dv/dz = 1.4 x 10 5 n(CO)/n(H 2 ) kms -1 pc _1 , respec- 
tively. These param eters are very clo se to those ob- 
tained previously by IWeifi et "all (|2007l ) in their single 
component model for the CO emission detected from this 
source. 

Adopting the same parameters for the water emit- 
ting region, we have computed the para-H 2 2 2 o — 
2n/CO J = 11 — 10 line flux ratio as a function 
of the assumed n(H 2 0)/n(CO) abundance ratio. In 
the case of H 2 0, the pumping of rotational transitions 
by far-infrared continuum radiation can strongly affect 
the predicted line fluxes. The importance of radia- 
tive pu mping in this source h as been discussed previ- 
ously by IRiechers et al.l (|2010|) for the case of HCN, al- 
though, in that case, pumping takes place through a 
low-lying vibrational band. Pumping through pure ro- 
tational transitions is relatively much more important 
for an asymmetric top molecule like water, because such 
molecules possess a more complex energy level structure 
than the simple ladder shown by spinless linear or di- 
atomic molecules (such as HCN and CO); furthermore, 
the lowest vibrational band of water lies at a consider- 
ably shorter wavelength (~ 6.3 /im) than that of HCN 
(~ 14.7 ^m), where the continuum radiation is consid- 
erably weaker. The dominance of radiative pumping in 
rot ational transitions of w a ter va por was also discussed 
by I Gonzalez- Alfonso et ail (|2010l ). in their recent analy- 
sis of the water line emission observed by Herschel toward 
the starburst galaxy Mrk 231. Under conditions where 
radiative pumping is dominant, the water line fluxes are 
almost inde pendent of the g as tem perature and density. 

Following IRiechers et al.l (|2010D . we define IRff as the 
sky covering factor of the infrared continuum source at 
the location of the molecular emission region; the mean 
intensity is then given by modified blackbody of the form 
IR ff (l -e- r )B l/ (220K), where B V {T) is the Planck func- 
tion and r = 1.05 (y/1 THz) 2 . In Figure [2] (upper left 
panel), we present the predicted para-H 2 2 2 n — 2n 
/ CO J = 11 — 10 line flux ratio, as a function of 
n(H 2 0)/n(CO) and for several different values of IRff: 
(red), 0.1 (magenta), 0.25 (blue), 0.5 (black), and 1.0 
(green). The results shown in Figure [5] clearly indicate 
the dominant role of radiative pumping; a detailed analy- 
sis indicates the importance of the In — > 2 2 q (rest frame 
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Fig. 2. — Molecular line ratios predicted by LVG models of APM 08279+5255, as described in the text. Labels in the lower two panels 
give the frequencies in GHz (wavelengths in fim) of water transitions that are predicted to be strong and potentially detectable. 



2.97 THz) transition in directly pumping the 2 2 o state of 
para-water, along with the 2o 2 — 33i (rest frame 6.45 THz) 
transition, which pumps the 331 state; the latter can de- 
cay subsequently to 220 • Dotted horizontal lines indicate 
the para-H 2 2 2 o — 2 n / CO J = 11 — 10 ratio measured 
in APM08279+5255 and its uncertainty. 

The upper right panel of Figure [5] shows entirely anal- 
ogous results for the lio — loi transition of ortho-water, 
with the horizont al dotted line indica ting the 3 a upper 
limit obtained bv iWagg et~aT1 (120061 ). A comparison of 
the results shown in the upper panels of Figure [5] indi- 
cates that a limited range of parameters is permitted by 
the measured value of the 2 2 o — 2n line flux and the 
upper limit on lio — loi- For the H 2 OPR of 3 as- 
sumed here, acceptable fits are obtained for IRff ~ 0.25 
and n(H 2 0)/n(CO) in the range - 0.003 - 0.03, al- 
though the range of acceptable parameters would obvi- 
ously broaden if OPR values smaller than 3 were permit- 
ted. In the regime of interest, the expected 2 20 — 2n 
line fluxes depend only weakly upon the H 2 abun- 
dance, the most important pumping transitions being 
optically-thick; thus, the exact range of acceptable values 
for n(H 2 0)/n(CO) depends strongly upon our estimate 
of the likely error in the measured line flux. Neverthe- 
less, the n(H 2 0)/n(CO) ratio inferred for an assumed 
OPR of 3 is apparently smaller than that typically mea- 
sured (~ 0.1 — 1) in hot core regions within our Galaxy 



(e.g. iBoonman fc van Dishoeckl (|2003l )). Observations of 
additional transitions will be needed to constrain the wa- 
ter OPR and abundance better. The lower panels of Fig- 
ure [2] present results for several other transitions that are 
potentially detectable from ground-based observatories, 
some of which are more strongly dependent upon the wa- 
ter abundance. The results shown in these panels were 
all obtained for IRg = 0.25 and OPR=3, and the labels 
indicate the rest frequencies and - in parentheses - red- 
shifted frequencies in GHz. Our results for transitions 
of para- and ortho-water appear, respectively, in the left 
and right panels. 

We note that our LVG solution for the H 2 Q emission in 
APM 08279+5525 not only matches that of iWeifi et all 
(|2007ft , using independent data, but the deduced velocity 
gradient dv / dz for the H 2 O emitting region is close to the 
expected vir i al valu e, as defined for example by eq. (5) of 
IGreve et all (|2009f ). K vir ~A (i.e. the dense gas emitting 
in the H 2 line emission is near virial equillibrium) . 

5. DISCUSSION 

The absence of detectable HF J = 1 — absorption 
in APM 08279+5255 is unexpected, given the low col- 
umn density of HF required to produce measurable ab- 
sorption. An important caveat in this analysis is the 
assumption of a Galactic HF/H 2 ratio that, in turn, is 
related to the elemental abundance of fluorine. Fluo- 
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rine nucleosynthesis — and thus the evolution of the flu- 
orine abundance in cosmic time — remain s poorly under- 
stood , with production in AGB stars fe.g. iCristallo et al.l 
120091 ) . in Wolf-Rayet stars (jMevnet fc Arnouldl2000t ) and 
neutrino- induced nucleosy nthesis in Type II supernovae 
(Woosley & Haxton 1998) all proposed as possible mech- 
anisms. While the face-on geometry of APM 08279+5255 
is not favorable for absorption studies, the lack of HF 
absorption is still puzzling and may be indicative of a 
lower fluorine abundance in this source compared with 
the Galactic ISM. Nevertheless, HF absorption may still 
prove to be a good tracer of H2 in high-redshift sources 
and additional observations of objects with different ge- 
ometries, over a wide range of redshifts, are urgently 
needed. 

Our LVG models indicate that the para-H 2 220 — 2n 
transition in AP M 8279+5255 is pre dominantly radia- 
tively pumped. IQmont et al.l (|2011l ) has reached sim- 
ilar conclusions regarding water excitation in J090302- 
014127B (SDP.17b) at z = 2.30. The para-H 2 2 20 -2 U 
line intensity in APM 8279+5525 is sensitive to the de- 
tails of the excitation model. Consequently, observations 
of this single transition do not provide a good estimate 
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of the water abundance. However, our LVG models sug- 
gests that many additional water lines should be de- 
tectable with the current millimeter-wave facilities. The 
transitions that are expected to be the strongest (see 
Fig. 2) are: 3 2 i — 3i 2 (rest frame frequency 1162.2 GHz), 
4 22 - 4i 3 (1207.6 GHz), 2 n - 2 02 ( 75 2.0 GHz), as well as 
two very high-energy transitions: 6 2 4 — 615 (1794.8 GHz), 
and 5 2 3 — 5i4 (1410.6 GHz), which are sensitive to the 
gas density. With multi-line observations, the excitation 
conditions and the water abundance will be much bet- 
ter constrained. This excitation scenario can further be 
tested with observations of the 2.97 THz pumping tran- 
sition, which is expected to appear in absorption. 
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